Quasinormal-mode spectrum of Kerr black holes and its geometric interpretation 
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There is a well-known, intuitive geometric correspondence between high-frequency quasinormal modes of 
Schwarzschild black holes and null geodesies that reside on the light-ring (often called spherical photon orbits): 
the real part of the mode's frequency relates to the geodesic's orbital frequency, and the imaginary part of the 
frequency corresponds to the Lyapunov exponent of the orbit. For slowly rotating black holes, the quasinormal- 
mode's real frequency is a linear combination of a the orbit's precessional and orbital frequencies, but the 
correspondence is otherwise unchanged. In this paper, we find a relationship between the quasinormal-mode 
frequencies of Kerr black holes of arbitrary (astrophysical) spins and general spherical photon orbits, which is 
analogous to the relationship for slowly rotating holes. To derive this result, we first use the WKB approxi- 
mation to compute accurate algebraic expressions for large-/ quasinormal-mode frequencies. Comparing our 
WKB calculation to the leading-order, geometric-optics approximation to scalar-wave propagation in the Kerr 
spacetime, we then draw a correspondence between the real parts of the parameters of a quasinormal mode 
and the conserved quantities of spherical photon orbits. At next-to-leading order in this comparison, we relate 
the imaginary parts of the quasinormal-mode parameters to coefficients that modify the amplitude of the scalar 
wave. With this correspondence, we find a geometric interpretation to two features of the quasinormal-mode 
spectrum of Kerr black holes: First, for Kerr holes rotating near the maximal rate, a large number of modes have 
nearly zero damping; we connect this characteristic to the fact that a large number of spherical photon orbits 
approach the horizon in this limit. Second, for black holes of any spins, the frequencies of specific sets of modes 
are degenerate; we find that this feature arises when the spherical photon orbits corresponding to these modes 
form closed (as opposed to ergodically winding) curves. 

PACS numbers: 04.25.-g, 04.30.Nk,4.70.Bw 



I. INTRODUCTION 

Quasinormal modes (QNMs) of black-hole spacetimes are 
the characteristic modes of linear perturbations of black holes 
that satisfy an outgoing boundary condition at infinity and an 
ingoing boundary condition at the horizon (they are the natu- 
ral, resonant modes of black-hole perturbations). These oscil- 
latory and decaying modes are represented by complex char- 
acteristic frequencies co — cor — icoj, which are typically in- 
dexed by three numbers, «, /, and m. The overtone number n 
is proportional to the decay rate of the perturbation, and / and 
m are multipolar indexes of the angular eigenfunctions of the 
QNM. 



A. Overview of quasinormal modes and their geometric 
interpretation 

Since their discovery, numerically, in the scattering of grav- 
itational waves in the Schwarzschild spacetime by Vishvesh- 
wara [1], QNMs have been thoroughly studied in a wide range 
of spacetimes, and they have found many applications. There 
are several reviews [2-6] that summarize the many discover- 
ies about QNMs. They describe how QNMs are defined, the 
many methods used to calculate QNMs (e.g., estimating them 
from time-domain solutions [7], using shooting methods in 
frequency-domain calculations [8], approximating them with 
inverse-potential approaches [9] and WKB methods [10, 11], 
and numerically solving for them with continued-fraction 
techniques [12, 13]), and the ways to quantify the excitation 
of QNMs (see, e.g., [14, 15]). They also discuss the prospects 



for detecting them in gravitational waves using interferometric 
gravitational-wave detectors, such as LIGO [16] and VIRGO 
[17], and for inferring astrophysical information from them 
(see, e.g., [18, 19] for finding the mass and spin of black 
holes using QNMs, [20, 21] for quantifying the excitation of 
QNMs in numerical-relativity simulations binary-black-hole 
mergers, and [22, 23] for testing the no-hair theorem with 
QNMs). There have also been several other recent applica- 
tions of QNMs. For example, Zimmerman and Chen [24] 
(based on work by Mino and Brink [25]) study extensions to 
the usual spectrum of modes generated in generic ringdowns. 
Dolan and Ottewill use eikonal methods to approximate the 
modal wave function, and they use these functions to study 
the Green's function and to help understand wave propagation 
in the Schwarzschild spacetime [26-28]. 

Although QNMs are well-understood and can be calculated 
quite precisely, it remains useful to develop intuitive and ana- 
lytical descriptions of these modes. An important calculation 
in this direction was performed by FeiTari and Mashhoon [9], 
who showed that for a Schwarzschild black hole, the QNM 
frequency (which depends only on a multipolar index / and an 
overtone index n) can be written in the eikonal (or geometric- 
optics) limit, I ^ 1, as 

CO ^ {I + l/2)Q.-iYL{n+ 1/2). (1.1) 

The quantities Q. and Yl are, respectively, the Keplerian fre- 
quency of the circular photon orbit and the Lyapunov expo- 
nent of the orbit, the latter of which characterizes how quickly 
a congruence of null geodesies on the circular photon orbit 
increases its cross section under infinitesimal radial perturba- 
tions [28, 29]. Equation (1.1) hints at an intriguing physi- 
cal description of QNMs: for modes with wavelengths much 
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shorter than the background curvature, the mode behaves as if 
it were sourced by a perturbation that orbits on and diffuses 
away from the light ring on the time scale of the Lyapunov 
exponent. 

Ferrari and Mashhoon [9] also derived an analogous result 
to Eq. (1.1) for slowly rotating black holes. They showed for 
If^m^ 1, the real part of the frequency is given by 

m 

i2 W a)orb + ^-qJYy^COprec, (1-2) 

where ©oit is now the Keplerian orbital frequency for the cir- 
cular photon orbit and COprec is the Lense-Thiring-precession 
frequency of the orbit (which arises because of the slow rota- 
tion of the black hole). The term proportional to (Bprec also has 
a simple geometric-optics interpretation. Inertial frames near 
the high-frequency wave at the light ring are dragged with 
respect to inertial frames at infinity, and this frame dragging 
causes the perturbation's orbit to precess about the spin axis of 
the black hole with a frequency, Oprec- If the orbit is inclined at 
an angle of sin^ 9 = /l{l + 1) (the ratio of angular momenta 
L? /L^ for quantized waves in flat space), then the projection 
of the precessional velocity onto the orbital plane scales the 
precessional frequency by a factor of ^ ;«/(/ + 1/2). 

Why the QNM frequency is multiplied by (/ + 1/2) is a 
feature that we will explain in greater detail in this paper In- 
tuitively, this term arises because the in the high-frequency 
limit, any wavefront traveling on null orbits will have an inte- 
gral number of oscillations in the 9 and (j) directions. For the 
wave to be periodic and single-valued, there must be m oscil- 
lations in the (j) direction. For the 9 direction, it is a Bohr- 
Sommerfeld quantization condition that requires / — \m\ + 1/2 
oscillations in this direction, which implies that there should 
be a net spatial frequency of roughly (/ + 1 /2). This increases 
the frequency of the radiation seen far from the hole by the 
same factor 

From this intuitive argument, we expect that the real part of 
the mode should be 

(Or^L [(Oorh + Y f^rec j , (1-3) 

where we define L = Z + 1/2. In this paper, we will show that 
an equation of the form of Eq. (1.3) does, in fact, describe the 
QNM frequencies of Kerr black holes of arbitrary astrophys- 
ical spins (and it recovers the result of Ferrari and Mashhoon 
for slowly spinning black holes). As we mention in the next 
part of this section, the exact details of the correspondence be- 
tween QNMs and photon orbits is richer for rapidly rotating 
black holes than for slowly rotating or static black holes. 

B. Methods and results of this article 

To derive Eq. (1.3) requires that we develop a geometric- 
optics interpretation of the QNMs of Kerr black holes with 
arbitrary astrophysical spins. Finding the correspondence re- 
quires two steps: first, we need to calculate the approximate 
frequencies using the WKB method; next, we must articulate 
a connection between the mathematics of waves propagating 



in the Kerr spacetime in the geometric-optics approximation 
and those of the WKB approximation (the first step). Finally, 
with the geometric-optics description of QNMs, we can make 
a physical interpretation of the spectrum (for example, the de- 
generacy or the lack of damping in the extremal limit). 

In Sec. II, we describe how we solve the eigenvalue prob- 
lem that arises from separating the Teukolsky equation [30] (a 
linear partial differential equation that describes the evolution 
of scalar, vector, and gravitational perturbations of Kerr black 
holes) into two nontrivial linear ordinary differential equa- 
tions. The two differential equations, the radial and angular 
Teukolsky equations, share two unknown constants — the fre- 
quency, 0), and the angular separation constant, A/,„ — that are 
fixed by the boundary conditions that the ordinary differen- 
tial equations must satisfy (ingoing at the horizon and outgo- 
ing at infinity for the radial equation, and well-behaved at the 
poles for the angular equation). The goal of the WKB method 
is to compute the frequency and separation constant approxi- 
mately. 

Although there has been work by Kokkotas [31] and Iyer 
and Seidel [32] using WKB methods to compute QNM fre- 
quencies of rotating black holes, their results were limited to 
slowly rotating black holes, because they performed an ex- 
pansion of the angular separation constant, A/,„, for small, di- 
mensionless spin parameters, a /M, and only applied the WKB 
method to the radial Teukolsky equation to solve for the fre- 
quency. In a different approach, Dolan developed a matched- 
expansion formalism for Kerr black holes of arbitrary spins 
that can be applied to compute the frequency of QNMs, but 
only for modes with / = \m\ and m = Q [27]. 

Therefore, it remains an outstanding problem to compute a 
WKB approximation to the quasinormal modes of Kerr black 
holes of large spins and for any multipolar index m. In Sec. II, 
we solve the joint eigenvalue problem of the radial and angular 
Teukolsky equations by applying a change of variables to the 
angular equation that brings it into the form of a bound-state 
problem in quantum mechanics. Applying the WKB method 
to the angular equation, we arrive at a Bohr-Sommerfeld con- 
dition that constrains the angular constant in terms of the fre- 
quency (and the indices / and m). Simultaneously, we can an- 
alyze the radial equation in the WKB approximation, and the 
two equations together define an system of integral equations, 
which can be solved for the eigenvalues. When we expand 
the Bohr-Sommerfeld condition in a Taylor series in terms of 
the numerically small parameter, a(o/l, the system of integral 
equations reduces to an algebraic system (which, in turn, leads 
to a simpler expression for the frequency). The approximate 
frequency agrees very well with the result that includes all 
powers of aco/l, and, in the eikonal limit, it is accurate to or- 
der 1 / 1 for Kerr black holes of arbitrary spins, for modes with 
any value of m, and for both the real and the imaginary parts 
of the frequency. 

To interpret the WKB calculation of Sec. II in the lan- 
guage of propagating waves in the geometric-optics limit 
within the Kerr spacetime, we analyze waves around a Kerr 
black hole in Sec. Ill using the geometric-optics approxima- 
tion and the Hamilton-Jacobi formalism. We confirm that 
the leading-order pieces of the WKB frequencies and angu- 
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Figure 1 : Low-overtone QNM spectrum of three Kerr black holes of different spins with approximate degeneracies in their spectra. From left to 
right, we plot the three lowest-overtone QNM excitations for (i) a/M = 0.69 in which {l,m) = {j,2) are black symbols and {I' ,m') = (_/ + 1, — 2) 
are blue symbols, where 7 = 3, ... ,9; (ii) a/M = 0.47 in which (l,m) = {j,3) are magenta symbols and {l',in') = (j+l,—3) are cyan symbols, 
where 7 = 3, ... ,9; (iii) a/M = 0.35 in which {l,m) = (y,4) are red symbols and {I' ,m') = {j+ 1,-4) are purple symbols, where 7 = 5, . . . , 10. 
For these spin parameters, the mode with positive values of m and (Or (a corotating mode) of index I is approximately degenerate with the 
mode with m' = — m, and cor (a counterrotating mode) of index /' = / + 1. 



lar constants correspond to the conserved quantities of the 
leading-order, geodesic behavior of the geometric-optics ap- 
proximation (specifically, the real part of O), the index m, and 
the real part of A/,„ are equivalent to the energy S', the z- 
component of the specific angular momentum L,, and Carter's 
constant i2 plus L^, respectively). The specific geodesies cor- 
responding to a QNM are, in fact, spherical photon orbits. The 
next-to-leading-order WKB quantities (the imaginary parts of 
(0 and A/„j) correspond to dispersive, wavelike corrections to 
the geodesic motion (they are the Lyapunov exponent and the 
product of this exponent with the change in Carter's constant 
with respect to the energy). Table I in Sec. Ill summarizes this 
geometric -optics correspondence. 

In Sec. IV, we make several observations about features of 
the QNM spectrum of Kerr black holes that have simple ge- 
ometric interpretations. First, we find that for extremal Kerr 
black holes, a significant fraction of the QNMs have a real 
frequency proportional to the angular frequency of the hori- 
zon and a decay rate that rapidly falls to zero; we explain this 
in terms of a large number of spherical photon orbits that col- 
lect on the horizon for extremal Kerr holes. Second, we ex- 
pand the WKB expression for the real part of the frequency as 
in Eq. (1.3), and we interpret these terms as an orbital and a 
precessional frequency of the corresponding spherical photon 
orbit. These two frequencies depend on the spin of the black 
hole and the value of m/L very weakly for slowly-rotating 
black holes, though quite strongly when the spin of the black 
hole is nearly extremal. Finally, we use the geometric -optics 
interpretation given by Eq. (1.3) to explain a degeneracy in 
the QNM spectrum of Kerr black holes, in the eikonal limit, 
which also manifests itself, approximately, for small / (see 
Fig. 1). The degeneracy occurs when the orbital and pre- 
cession frequencies, cOoit and CDprec are rationally related (i.e., 
('^oib I Oprec — p / q fo^ integers p and q) for a hole of a specific 
spin parameter, and when the corresponding spherical photon 
orbits close. By substituting this result into Eq. (1.3) one can 
easily see that modes with multipolar indexes I and m become 
degenerate with those of indexes I' — I + kq and m' — m — kp 
for any non-negative integer k, in the eikonal limit (note that 
in Fig. 1 , we show an approximate degeneracy for ^ = I and 
for three spin parameters, such that q/ p = 1/4, I /6, and 1/8, 



respectively.) 

To conclude this introduction, we briefly summarize the or- 
ganization of this paper: In Sec. II, we review the Teukol- 
sky equations, and we then describe the WKB formalism 
that we use to calculate an accurate approximation to the an- 
gular eigenvalues A/,„ = Af^^^ + iA\^^ and QNM frequencies 
(0 = (Or — i(Oj, in the eikonal limit L ^ 1 and for holes of 
arbitrary spins. We verify the accuracy of our expressions in 
Sec. II D by comparing the WKB frequencies to exact, numer- 
ically calculated frequencies. In Sec. Ill, we develop a corre- 
spondence between the WKB calculation and mathematics of 
wave propagation within the Kerr spacetime in the geometric- 
optics limit, using the geometric-optics approximation and 
the Hamilton-Jacobi formalism. At leading-order, the QNM 
frequencies and angular eigenvalues correspond to the con- 
served quantities of motion in the Kerr spacetime for spher- 
ical photon orbits; at next-to-leading order in the geometric- 
optics approximation, we connect the the decaying behavior 
of the QNMs to dispersive behaviors of the waves. Finally, in 
Sec. IV, we interpret aspects of the QNM spectrum geomet- 
rically, such as the vanishing of the damping rate for many 
modes of extremal black holes, the decomposition of the fre- 
quency into orbital and precessional parts, and the degenera- 
cies in the QNM frequency spectrum. Finally, in Sec. V, we 
conclude. We use geometrized units in which G = c = 1 and 
the Einstein summation convention throughout this paper 



II. WKB APPROXIMATION FOR THE 
QUASINORMAL-MODE SPECTRUM OF KERR BLACK 
HOLES 

In this section, we will derive expressions for the frequen- 
cies of quasinormal modes of Kerr black holes using the WKB 
approximation. We will need to compute the real and imagi- 
nary parts to an accuracy of (9(1) in terms of Z 1, which im- 
plies that we must calculate (Or to leading and next-to-leading 
order and (Oj to leading order Here, we will focus on obtain- 
ing an analytic approximation to the frequency spectrum, and 
we will leave the geometrical interpretation of our results until 
the next section. 
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Before specializing our results to the angular and radial 
Teukolsky equations, we will review a basic result about the 
WKB expansion that we will use frequently throughout this 
paper; a more complete discussion of WKB methods can be 
found in [11]. Given a wave equation for ^/(x) 



(2.1) 



we will expand the solution as i/a = g5'o/e+Si+£52+...^ where the 
leading and next-to-leading action variables are given by 



5o = ±/ / \JlJ{x)dx , 



5i = --logt/W. 



(2.2a) 
(2.2b) 



The formulas above will be the basis for our analysis of the 
radial and angular Teukolsky equations in the next sections. 



A. The Teukolsky equations 

Teukolsky showed that scalar, vector, and tensor perturba- 
tions of the Kerr spacetime all satisfy a single master equation 
for scalar variables of spin weight s\ moreover, the master 
equation can be solved by separation of variables [30]. We 
will use M to denote our scalar variable, and we will separate 
this scalar wave as 



M(f,r,0,0)=e-""V'"'''M,(r)Me(0). 



(2.3) 



B. The angular eigenvalue problem 

We will first find an expression for A/„, in terms of O), /, and 
m, by analyzing the angular equation in the WKB approxima- 
tion. By defining 



X = log tan 







(2.6) 



and dx — dQ I sin 0, we can write the angular equation as 



dx^ 



Vmq -0, 



(2.7a) 



where 



^fl^CO^cos^esin^e-mHAtasin^e. (2.7b) 

When written in this form, it is clear that uq must satisfy a 
boundary condition that it be as jc — >^ ±oo (which corresponds 
to ^ 0, k). Furthermore, the angular equation is now in a 
form that is amenable to a WKB analysis (which will be the 
subject of the next part). 

First, however, we outline how we will perform the calcu- 
lation. Because the frequency O) = Or — /O/ is complex, the 
angular eigenvalue A/,,,, a function of o, must also be com- 
plex. We will write 



A;,„ =A 



Im 



iA 



(2.8) 



to indicate the split between real and imaginary parts. We will 
treat a real-valued CO — (Or in the first part of this section, and, 
therefore, a real-valued A f^^{cOR); we shall account for — /©/ 
by including it as an additional perturbation in the next part of 
this section. 



Then, at the relevant order in / ^ 1 , the angular equation for 
ug{0) can be written as 



1. Real part ofAi,„ for a real-valued CO 



1 



sine dO 









+ 







fl^OJ^cos" 9 - 



sin^e 



lie =0, 
(2.4) 



where A/„, is the angular eigenvalue of this equation. The 
equation obeyed by the radial function Ur{r) is 



drl + (r2+a2)2«'--^' 



d 



+ dr 



(2.5a) 



with 



K 



A/^„ = A/m + fl^O)^ — lamco , 
A = r2-2Mr + fl2. 



(2.5b) 
(2.5c) 
(2.5d) 

Here we have used the facts that (Or ^ 0(1), (Oj ^ 0{l), 
m ~ 0{l) to drop terms that are of higher orders in the expan- 
sion than those that we treat. Note that the spin s of the per- 
turbation no longer enters into these equations after neglecting 
the higher-order terms. 



For (Or e M, we will compute the eigenvalues Af^^^{cOR), of 
Eq. (2.7a) for standing-wave solutions that satisfy physical 
boundary conditions. At the boundary, 9 = Q,n (or x = 
the potential satisfies — 



4^ 



~m independent of the value of 
this implies that the solutions to Eq. (2.7a) behave like 
decaying exponential functions at these points (i.e., the wave 
does not propagate). For there to be a region where the so- 
lutions oscillate (i.e., where the wave would propagate), A/,„ 
must be sufficiently large to make > in some region. De- 
pending on the relative amplitudes of A/,„ and a^(0^, either 
has one maximum at = 7i/2 (when A/,„ > a^(0^), or two 
identical maxima at two locations at symmetrically situated 
around 9 = n/2 (when A/,„ < a^(0^). It turns out that the re- 
gion where the maximum of > is centered around n/2; 
therefore, all solutions fall into the former category rather than 
the latter. 

The length scale over which the function ug varies is 
1 /W^, and the WKB approximation is valid only if the po- 
tential does not vary much at this scale. Therefore, to use 
the WKB approximation, we require that 



1 dV*" 



/ye d9 



(2.9) 
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This condition applies regardless of the sign of . Empiri- 
cally, we find this condition to hold for in Eq. (2.7a), ex- 
cept around points at which = 0. We will refer to these as 
turning points, and they can be found by solving for the zeros 
of the potential. 



a^a)| cos^ 0±sin^ 0± 



which are given by 



^Af„sin2 0±^O, (2.10) 



sin-^ 0± = 



2m^ 



A,,„ + fl2co2 ip / (A,,„ + fl2c(j2 )2 +4^2 



:. (2.11) 



Using the leading and next-to-leading WKB approximation, 
we can write the solution to the wave equation in the propaga- 
tive region, x_ < x < x^, as 



Mg(x) = — , (2.12) 

where a± are constants that must be fixed by the boundary 
conditions that the solution approach zero at = 0, K. For 
X > x+, we find 



ue{x) = 



1/4 



(2.13a) 



and X < 



ue{x) 



-J.^-dx'^-ve{x') 



[^(x) 



1/4 



(2.13b) 



with c± also being constants set by the boundary conditions. 
Note that outside of the turning points, we have only allowed 
the solution that decays towards x — > ±°o. 

Around the turning points x±, the WKB approximation 
breaks down, but ug can be solved separately by using the fact 
that Vg{x ^ x±) oc X — x±. Solutions obtained in these regions 
can be matched to Eqs. (2.12)-(2.13b); the matching condi- 
tion leads to the Bohr-Sommerfeld quantization condition 



then all three of these quantities are 0(1) in our expansion 
in L. From these definitions, we can re-express the limits of 
integration as 



sin^ e± = 

and the integral as 



(2.17) 



sin^e 



2^a^cos^e^{\-\li\)n. (2.18) 



For each set of quantities {aR,lJ.,Q.R), we can express Ur 
as an implicit function involving elliptic integrals; however, if 
we treat aQ.R as a small parameter, then the first two terms in 
the expansion are 



2o2 



afisa 1 



-(1- 



(2.19) 



We derive and discuss this approximation in greater detail in 
Appendix A. For a = 0, we note that this is accurate with a 
relative error of (9(1 /L^), because for a Schwarzschild black 
hole 



Af„fw = /(/ + l)-.(.+ l). 



(2.20) 



As we will confirm later in Figs. 2 and 4, Eq. (2.19) is an 
excellent approximation even for highly spinning black holes. 

To understand intuitively why the approximation works so 
well, we will focus on corotating modes (i.e., those with posi- 
tive and large m, or ji near unity), which have the highest fre- 
quencies and, therefore, the largest possible values for oCIr. 
For a fixed value of {l,m), cOr is a monotonically increasing 
function of a, and 



a)^'"(fl) < <'(fl =M)= mQ."ij 



a=l 



2M 



(2.21) 



In setting this upper bound, we have used the result that the 
low-overtone QNM frequencies approach mQ.H for m > and 
for extremal black holes (first discussed by Detweiler [33]); 
we have also used CIh to denote the horizon frequency of the 
KeiT black hole. 



d9\ a^col cos2 0- 



sin^e 



Al^^{L~\m\)K. (2.14) 



Here we have defined 



L=l 



(2.15) 



which will be used frequently throughout this paper. The lim- 
its of the integration are the values of Q where the integrand 
vanishes [the turning points of Eq. (2.11)]. 
If we define 



;U=-, aR{a,iJ.) = ^, Q.R{a,n) = — 



(2.16) 



2Mr4 



(2.22) 



and r+ to indicate the position of the horizon [note that r+ {a - 
M) ~ M]. Normalizing Eq. 2.21 by L, we find 



aQ.R < {ll/2){a/M) < 1/2. 



(2.23) 



Even for the upper bound uSIr = 1 /2, the relative accuracy of 
Eq. (2.19) is still better than 0.2%. 



2. Complex Aj^ for a complex CO 

To solve for the next- to-leading-order corrections to A/,„, 
we must compute the imaginary part A^^^^. Because coj ^ cor. 



when we allow co = (Or — icoj to be a complex number in the 
angular eigenvalue problem (2.4), we can treat the term linear 
in coj as a perturbation to the angular equation. Using the 
perturbation theory of eigenvalue equations, we find that 



A|,„ = -2a^ cOrCOj {cos^ 6) , 



(2.24) 



where 



(cos^ 9) 



cos^ 9\ug\^ sinddO 



\ue\ sinOdO 

cos^ 



--dd 



Im 



(2.25) 



-de 



fl2a)|cos2 9 — 



By taking the derivative of both sides of the Bohr-Sommerfeld 
condition (2.14) with respect to the variable z = oCOr and by 
treating A/,„ as a function of z, we can rewrite the above ex- 
pression as 



(cos^e) 



1 dAliz) 



Z=C10)R 



2z dz 

Substituting this expectation value into Eq. (2.24), we find 



(2.26) 



= ami 



dz 



(2.27) 



Equation (2.27) defines a numerical prescription for comput- 
ing A/„, = Af^^ + iA\^^. This approach is quite natural: as O) 
becomes complex, A/^ is the analytic function whose value 
on the real axis is given by Af^^. The approximate formula 
(2.19), therefore, becomes 



L2 



(2.28a) 



1 /I I |2^ 



if we define 

V r, (O) = (2^ 2M 

(2.29b) 

Note here that V is an analytic function of o, and that it is 
real-valued when O) is real. 

In general, the WKB approximant for Ur is given at leading 
order by 

= b+e'^^'' V^)dr', + b_e-'r* V^)dri ^ (2.30) 

although in order to obtain a mode which is outgoing at ^• 
+00 (the same as r oo) and ingoing at —°° (r r+), 
we must have 



for the region containing r +°o, and 



(2.31a) 



(2.31b) 



for the region containing — > —0°. Intuitively speaking, a 
solution to Eq. (2.29a) will satisfy the asymptotic behavior 
above if V « around a point r = ro, and > on both sides. 
Then, the WKB expansion (2.30) is valid in the two regions 
on both sides of r = ro, and the solution in the vicinity of ro 
must be obtained separately by matching to the WKB approx- 
imation. The matching will constrain the frequency, thereby 
giving a method to determine (O. A detailed calculation of 
this procedure has been carried out by Iyer and Will [11] to 
high orders in the WKB approximation; the only difference 
between our calculation and their result at lower orders comes 
from the more complex dependence of V on (O in our case 
(particularly because A/,„ depends on O) in a more involved 
way). 



1. Computing (Or 

From Iyer and Will [11], the conditions at the leading and 
next-to-leading order that must be solved to find (Or are 



V''(ro,«R) 



dV 



0. 



(2.32) 



[ro.DR) 



(2.28b) After a short calculation, these conditions can be expressed as 



for a complex frequency co. 



C. The radial eigenvalue problem 

Now that we have solved for the angular eigenvalues A/,„ in 
terms of (O, we turn to the radial Teukolsky equation. From 
Eq. (2.5a), we see that the radial equation is already in the 
form 



d Ur 

~dri 



(2.29a) 



rg +0-^ + fl^ 



»4 



7W) 



where we have defined 

P{z) = Ja{z)+z^-2iiz 



l + --2^iz- 



(2.33a) 
(2.33b) 

(2.34a) 
(2.34b) 
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In deriving Eq. (2.33b), we have used the fact that at r > r+, 
(r^+fl^)^/A is a monotonically increasing function, and, 
therefore the extrema of V are the same as those of V (r^ + 
fl^)^/A; we then also used the fact that the quantity within the 
square brackets in Eq. (2.33b) is always nonzero at points at 
which V = 0. 

One method of jointly solving Eqs. (2.33a) and (2.33b) 
would be to use Eq. (2.33b) to express Q.r in terms of ro 

{M-ro)na 



(ro-3M)r2 + (ro+M)a2 



(2.35) 



and then insert this into Eq. (2.33a) to obtain ro; finally Q.r 
can be obtained by substituting this ro back into Eq. (2.35). 
If we use the approximate formula (2.34b) in this process, the 
equation for ro becomes a sixth-order polynomial in x = ro/M, 
the roots of which can be found efficiently. For conveninece, 
we write this polynomial here 



2x'*(x-3) 



-4x^[{l-li^)x^-2x-3{l-li^)]{a/Mf 



+ {l-^'-)[{2-^^)x^+2i2 + n^)x+{2^ 



-li^)]{a/Mf. 

(2.36) 



For each pair {ii,a/M), there are in general two real roots for 
X, which correspond to the two possible values of rg/M (and 
the two real frequencies with opposite signs). 

Note that the procedure above will not work when m = 
(when both the numerator and denominator of Eq. (2.35) van- 
ish). In this case, we can directly require that 

{r„-3M)rl + {rp+M)a^ =0. (2.37) 

The solution, rp can be inserted into Eq. (2.33a) and the result 
can be expressed in terms of elliptic integrals 



i2fl(fl,M-0)-±- 



(r2+fl2)EllipE 



a^A(rp) 
2\2 



-a2) 



(2.38) 



where EUipE denotes an elhptic integral of the second kind. 
Here we have used the subscript p for this special case, be- 
cause this mode will turn out to correspond to polar orbits. 

We plot in Fig. 2 the relative error in Q.r that comes from 
using the approximate expression for A/„, [Eq. (2.28a)] rather 
than the exact Bohr-Sommerfeld condition. The error is al- 
ways less than ^ 10^"* (we scale the quantity plotted on the 
vertical axis by 10^), and therefore, we will use the approx- 
imate expression for A/,„ for the remaining calculations in- 
volving Q.R throughout this paper In Fig. 3, we plot D.r for 
fl/M = 0, 0.3, 0.5, 0.9, 0.99, and 1 (the flat curve corresponds 
to a/M — 0, and those with increasing slopes are the increas- 
ing values of a/M). While for low values of a/M below ^ 0.5, 
Q.R depends roughly linearly upon /i, for higher values of a/M 
(and for pL > 0), CIr grows more rapidly than linearly. For 
a/M = 1, Q.R — > 1/2 when /i — > 1, as anticipated. 



2. Computing (Oj 

At leading order, the imaginary part (Oj can be calculated 
using the same procedure set forth by Iyer and Will [11]. They 
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Figure 2: Difference in fl«(a,/i) [Eq. (2.35)] that arises from us- 
ing the approximate formula for A/,„ [Eq. (2.28a)] as opposed to the 
exact formula. Here a/M = 0.7, 0.9,0.95, and 0.99 correspond to 
black solid, red dashed, blue dotted, and purple long-dashed curves, 
respectively. The quantity plotted on the vertical axis has been scaled 
by 10^ 
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Figure 3: Real part of the QNM spectra from the WKB approxi- 
mation. Black solid curves show fl for a/M = (the flat curve) 
and a/M = 1 (the curve that increases towards 0.5); red (light gray) 
dashed and dotted curves show a/M = 0.3 and 0.5, while blue (dark 
gray) dotted and dashed curves show a/M = 0.9 and 0.99. 



find that 



CO/ = -(n-f 1/2) 



(«+l/2)n/(a,^) 



(2.39) 
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In our calculation, we must also take into account that V also on O). If we use the approximate formula for a, we obtain a 
depends on co through the angular eigenvalue's dependence reasonably compact expression for 



(2.40) 
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Figure 4: Difference in £lj{a,ix) [Eq. (2.39)] from using the approx- 
imate formula for A;„, [Eq. (2.28a)] rather than the exact formula. 
Here a/M = 0.7, 0.9, 0.95, and 0.99 correspond to black solid, red 
dashed, blue dotted, and purple long-dashed curves, respectively. We 
scale the quantity plotted along the vertical axis by 10"^ in this fi sure. 



In Fig. 4, we plot the relative error in i2/ from using the 
approximate formula for A/,„ identically to that in Fig. 2 (al- 
though here we scale the quantity plotted on the vertical axis 
by 10'*). Because the error is always less than ^ 10^^, we will 
use the approximate expression for A/„, for computing Q.j in 
the remainder of this paper. In Fig. 5, we plot i2/(fl,/x) for 
several values of a/M (the same as those in Fig. 3, though not 
a/M = 0.3). The curve for a/M = is flat, and those with 
larger spins have more rapidly decreasing slopes for increas- 
ing values of a/M. It is interesting to note that in the limit 
fl — ^ 1, Q.I becomes zero for values of /i in a finite interval 
0.74 < jU < 1 (not only for = 1 does Q.j vanish). We will 
put forward an explanation for this phenomenon in Sec. IV, 
after we make connections between QNMs and wave propa- 
gation in the Kerr space time. 



D. Accuracy of the WKB approximation 

Because we calculated the leading and next-to-leading or- 
ders in the WKB approximation to cor, we expect that the rel- 
ative error for increasing L should scale as 0{1/L^). For the 
imaginary part, we computed only the leading-order expres- 
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Figure 5: Imaginary part of the QNM spectrum computed in the 
WKB approximation. Black solid curves show Q.j for a/M = 
(again the flat curve) and a/M = 1 the curve that decreases and heads 
to zero. The red dashed curve shows a/M = 0.5, while blue dot- 
ted and dashed curves show a/M = 0.9 and 0.99, respectively. For 
a/M = 1, modes with jl > 0.74 approach zero (modes that do no 
decay), while others still decay. 



sion, and we would expect that the relative eiTor might scale 
as 0{l/L). In addition, because at this order of approxima- 
tion, we do not account for the spin of the wave, we anticipate 
that the error for the gravitational modes may be larger than 
those for scalar modes. In Figs. 6-9, we confirm most of these 
expectations, but we find the somewhat unexpected result that 
the relative eiTor for the imaginary part also scales as 0{1/L^). 

In Fig. 6, we compare the WKB approximation to cor 
with numerical computations of the i = 2, gravitational-wave, 
quasinormal-mode spectra; specifically, we plot the fractional 
error against = m/L, for Z = 2, 3, . . . , 14, and for black holes 
of spins a/M = 0.3, 0.5, 0.9, and 0.95. The relative error 
clearly converges to 0{L^). Even for / = 2, the relative error 
tends to be < 30%, and at Z > 3 the relative error stays below 
^ 1.5L^2 (this means error is < 10% for / = 3 and higher). 

In Fig. 7, we compare the WKB spectra with i = scalar 
quasinormal-mode spectra, for the same values of / and the 
same black-hole spins. We find a much better agreement. For 
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Figure 6: Fractional error, Scor/cor, of the WKB approximation to the s = 2, gravitational-wave, quasinormal-mode spectrum, multiplied by 
L?. The four panels correspond to four different spins which (going clockwise from the top left) are a/M = 0.3, 0.5, 0.95, and 0.9. Errors for 
/ = 2,3,4 are highlighted as red solid, brown dashed, and pink dotted lines, while the rest (/ = 5, . . . , 14) are shown in gray. This shows that 
the relative error approaches the 0{l/L?) scaling quite quickly. 
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Figure 7: Fractional error, ScO/t/cog, of the WKB approximation to the i- = 0, scalar- wave, quasinormal-mode spectrum, again scaled by L?. 
The four panels correspond to the same four spins in Fig. 6. The points shown in the four panels are for values of I in the range / = 2, 3, . . . , 14. 
Because all values of / nearly lie on the same curve, the relative error has converged at an order 0{l/L?) even for very low /. The overall error 
is also significantly lower than that for the i- = 2 modes. 
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Figure 9: Fractional error, 5cOj/cOj, of the WKB approximation to the i = 0, scalar-wave, quasinormal-mode spectrum, again multiplied by 
L?. The four panels and the points are shown in the same way as in Fig. 7, and there is a similar rapid convergence of the error. 
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all Z > 2 modes, the relative error stays below 4 x lO^^L^^. 
This suggests that coupling between the spin of the wave (i.e., 
its tensor polarization) and the background curvature of the 
Kerr black hole is the main source of error in our WKB ap- 
proximation. 

In Figs. 8 and 9, we perform the same comparisons as in 
Figs. 6 and 7 for the imaginary part of frequency. Surprisingly, 
we find that for both s = and 2, the relative error in coj is 
0{L^^). For s ~ 0, the relative error is < 6 x lO^^L"^, while 
for s — 2, the error is < L^^. 

With this comparison, we conclude our direct calculation of 
the QNM spectrum by WKB techniques. We will discuss ad- 
ditional features of the QNM spectrum in Sec. IV, but before 
doing so, we will develop a geometric interpretation of our 
WKB results. Doing so will help us to develop more intuition 
about our WKB expressions. 

III. GEOMETRIC OPTICS IN THE KERR SPACETIME 

In this section, we first briefly review the formalism of geo- 
metric optics, which describes the propagation of waves with 
reduced wavelengths X that are much shorter than the space- 
time radius of curvature, R, and the size of the phase front. 
In the geometric-optics approximation, the phase of the 
waves remains constant along null geodesies (rays), while the 
amplitude can be expressed in terms of the expansion and con- 
traction of the cross-sectional area of bundles of null rays. 
We will then specialize the geometric -optics formalism to the 
Kerr spacetime, and we will write down the most general form 
of propagating waves in the geometric-optics approximation. 
Using the Hamilton-Jacobi method, we see that the waves' 
motion can be related to the null geodesies in the spacetime. 
By applying boundary conditions to the approximate wave, 
we obtain expressions for the quasinormal-mode waveforms 
and their corresponding complex frequency spectra and angu- 
lar separation constants, in the eikonal limit. 

A. Geometric optics: general theory 

Here we briefly review the geometric-optics approximation 
to scalar-wave propagation (see, e.g.. Section 22.5 of Ref. [35] 
for details). A massless scalar wave u propagating in curved 
spacetime satisfies the wave equation 

gA'^V^Vv« = 0. (3.1) 

If we write 

M=Ae'*, (3.2) 
then at leading order in A / we have 

^^%fev = 0, k^ = d^^, (3.3) 
while at next-to-leading order, 

2^''5^1ogA + V^^^ =0. (3.4) 



Note that Eq. (3.3) also implies that is geodesic, 

'fe'^V^/tv =/t''Vj,Vv4> = /t''VvV^4> = fc''Vv/tj^ =0. (3.5) 

Equations (3.2)-(3.5) encode information about the trans- 
port of the amplitude A and phase <I> along a null geodesic (or 
a my). The phase should be kept constant, because Eq. (3.3) 
states 

k^d^^ = Q, (3.6) 

while the amplitude is transported along the ray in a manner 
that depends upon the propagation of neighboring rays. Be- 
cause the 2D area, jz/, of a small bundle of null rays around 
the central ray satisfies the equation 

y^k^^k^d^logs^, (3.7) 

it is possible to show from Eq. (3.4) that 

kf'dy.^s^^l^A^^O, (3.8) 

which implies A oc s^^^l'^. 

The transport equations (3.6) and (3.8) provide a way to 
construct a wave solution from a single ray; therefore, any 
solution to the wave equation (3.1) in a 4D spacetime re- 
gion can be found from a three-parameter family of null rays 
(with smoothly varying initial positions and wave vectors) by 
assigning smoothly varying initial values of (OjA) and then 
transporting these values along the rays. (We use the phrase 
"smoothly varying" to mean that the values of ('t'jA) must 
change on the scale of ^ ^ A .) We note it is often convenient 
to divide the three-parameter family of initial positions of the 
null rays into two-parameter families of rays with constant ini- 
tial values of 4>. The constant-<I> surfaces are the initial phase 
fronts, which, upon propagation along the rays, become 3D 
phase fronts of the globally defined wave. The more usual 2D 
phase fronts, at a given time, are obtained if we take a particu- 
lar time slicing of the spacetime and find the 2D cross sections 
of the 3D phase fronts in this slicing. 

The above formalism describes wave propagation up to 
next- to-leading order in ^/A, which will be enough for us to 
build a geometric correspondence for both the real frequency, 
the decay rate, and the angular separation constant of QNMs 
in the Kerr spacetime. 

B. Null geodesies in the Kerr spacetime 

Now let us review the description of null geodesies in the 
Kerr spacetime using the Hamilton-Jacobi formalism. In gen- 
eral, the Hamilton-Jacobi equation states 

^^'(<3m^)M-0, (3.9) 

where S{x^^ is called the principal function. For the Kerr 
spacetime, the Hamilton-Jacobi equation can be solved via 
separation of variables (see, e.g., [34]), through which the 
principal function can be expressed as 

5(f,0,0,r)=50(0)+L,0+5,-(r)-^f, (3.10) 
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where and are constants that are conserved because of the 
the timelike and axial Killing vectors of the Kerr spacetime. 
Physically, <§ and represent the energy and z-directed spe- 
cific angular momentum of the massless scalar particle. The 
functions and Se{d) are given by 



Seie) 



^@{e')de' , 

(3.11a) 



where ^%{r) and 0(0) are given by 

^(r) = [S{r^ +a^)-L,af -^[{L-,-aSf + ^], (3.11b) 

0(0) = ^-COs2 0(L2/sin2 0-fl2^2-)^ (1,11c) 

and A(r) is given in Eq. (2.5d). The constant ^ is the Carter 
constant of the trajectory, which is a third conserved quantity 
along geodesies in the Kerr spacetime. 

The principal function S{x^ ;<S',L^,J2) contains information 
about all null geodesies; equations of motion for individual 
null geodesies are given by first choosing a particular set of 
((fjL^, cS), and then imposing 



^.0 ^=0 ^=0 



(3.12) 



These conditions lead to at set of first-order differential equa- 
tions 



dt + c?- 



dX 
d^ 



[f?(r2 + fl2)-Lj.a] -a{aS'5i\\^e-L,) 



(3.13a) 



L, \ , a[S{r^ + a^)-L,a\ 

^aS-^^\+^ (3.13b) 

dX \ sm^ / A 



dr ^j-^^ dO 
dX ' dX 

where we have defined 



0, 



(3.13c) 



d ^ d ^ 9 9 7„ 
-^ = I,—r, E = r2 + fl2cos2 0, 
dX dQ 

and C, is an affine parameter along the null geodesies. 



(3.14) 



C. Correspondence with quasinormal modes 

Here we will find connection between the general set of 
wave solutions in the previous section, and the particular 
solutions that correspond to a quasinormal modes, in the 
geometric -optics limit. Specifically, we will look for waves 
that propagate outwards at infinity and down the horizon. 
With this correspondence, we will be able to make a geometric 
interpretation of our WKB results from Sec. 11. 



I. Leading order: conserved quantities of rays and the real parts 
of quasinormal-mode parameters 

It is straightforward to note that the Hamilton- Jacobi equa- 
tion (3.9) is identical to the leading-order geometric -optics 



equations, if we identify the phase, <I>, with the principal func- 
tion, S. Therefore, at leading order, we can write 

-i<gt JL-J +iSs +iS, 



JS . 



(3.15) 



where we recall that the amplitude A differs from unity only 
at next-to-leading order (we will treat it in the next subsec- 
tions). Here, we have a chosen set of conserved quantities, 
((f,^,L,), to identify the wave we wish to connect with a 
quasinormal-mode solution. The region in which the wave 
propagates is identical to the region in which geodesies with 
these conserved quantities can propagate. In addition, for each 
point in this region, there is one (and only one) geodesic pass- 
ing through it; that we have ± in front of Se and Sr means 
only that either propagation direction could be a solution to 
the wave equation. 

Now we note that u, a scalar wave in the Kerr spacetime, 
must separate as in Eq. (2.3). By comparing Eq. (2.3) and 
Eq. (3.15), we can immediately identify that 



(3.16) 



Because is a real quantity (the conserved energy of the null 
geodesic), we see that at leading order, the wave does not de- 
cay. Next, we also observe that in order for u to be consistently 
defined in the azimuthal direction, (of the null geodesies 
that S describes) must be an integer This allows us to make 
the second identification 



(3.17) 



Comparing Sq from Eq. (3.11a) and uq from Eqs. (2.12) 
and (2.7b) (focusing on one direction of propagation, and 
ignoring next-to-leading-order terms), we can also conclude 
that 



(3.18) 



At this stage, given any set of {S',^,L^), we will be able 
to find a wave solution that exists in the region in which the 
geodesies travel. Not all such sets of conserved quantities cor- 
respond to quasinormal modes, however, because they may 
not satisfy the correct boundary conditions of QNMs. 

We will first explain the conditions on the radial motion 
of the geodesies that will allow these particular geodesies to 
correspond to a wave that satisfies outgoing and downgoing 
conditions at ±°o, respectively. If the radial geodesies 
satisfy £i > everywhere, then there will be traveling waves 
across the entire axis, which will not satisfy the boundary 
conditions; if there are two disconnected regions of traveling 
waves, however, waves will scatter off the potential on each 
side, and they will also fail to satisfy the boundary conditions. 
The only way to satisfy the boundary conditions is to have 
a point ro at which ^ = and = 0, in which case there 
will be a family of geodesies on each side of r = ro (with 
each member a homoclinic orbit which has r — > ro on one end) 
and a spherical orbit with constant r = ro. The corresponding 
wave has zero radial spatial frequency at r = ro, and this fre- 
quency increases towards r < ro and decreases towards r > ro. 
Noting that 



V{r'+a'f 



(3.19) 
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Table I: Geometric-optics correspondence between the parameters of a quasinormal mode, (co, Ai,„, I, and m), and the conserved quantities 
along geodesies, L~, and To establish a correspondence with the next-to-leading-order, geometric-optics approximation, the geodesic 
quantities if and J2 must be complex. 



Wave Quantity 


Ray Quantity 


Interpretation 




S 


Wave frequency is same as energy of null ray 
(determined by spherical photon orbit). 




L, 


Azimuthal quantum number coiTesponds to z angular momentum 




(quantized to get standing wave in direction). 






Real part of angular eigenvalue related to Carter constant 
(quantized to get standing wave in 6 direction). 


COi 




Wave decay rate is proportional to Lyapunov exponent 




of rays neighboring the light sphere. 


A' 




Nonzero because £0/ 7^ 
(see Sees. II B 2 and III C 3 for further discussion). 



the condition 



^ = ^' = 



(3.20) 



is the same as the condition, Eq. (2.32), which determines (Or 
in terms of L and m in the WKB approximation. It is worth 
mentioning that although the condition of Eq. (3.20) imposed 
on {(oj^jL^) can be interpreted most easily as the condition 
for a spherical photon orbit, the wave function for the quasi- 
normal mode we are considering is not localized around that 
orbit. The wave function at leading order, in fact, has a con- 
stant magnitude at every location that homoclinic orbits reach 
(i.e., the entire r axis). We will derive the amplitude correc- 
tions in the next section. 

The quantization of the frequency cor in terms of the multi- 
polar indices I and m arises from the quantization of the mo- 
tion in the angular directions. For the azimuthal direction, it is 
easy to see that for the wave function to be single-valued, we 
need to impose = m eZ. For the 9 direction, we note that 



& = V' 



' sin^ . 



(3.21) 



and the 0-quantization condition for the wave, Eq. (2.14), 



IS 



@d9 = [L- \m\)%. 



(3.22) 



This corresponds to the Bohr-Sommerfeld condition for a par- 
ticle moving in a potential given by 0. Consequently, the con- 
dition for a standing wave along the 9 direction (at leading 
order) is equivalent to 



i? — Aii„ {(Oro) — 

9 a COJ 



m 

Z2 



.21 



(3.23) 



In summary, we connected the QNM's wave function to 
the Hamilton-Jacobi principal function of homoclinic null 
geodesies (at leading order). These geodesies have the same 
energy. Carter constant, and z-component of its angular mo- 
mentum as a spherical photon orbit; however only spherical 



orbits with quantized Carter constants and z-angular momenta 
correspond to quasinormal modes. In Table 1, we summarize 
our geometric -optics correspondence; so far we have identi- 
fied the first three entries on the table. We can find the next two 
correspondences by investigating next-to-leading-order geo- 
metric optics in the next part. 



2. Next-to-leading order: radial amplitude corrections and the 
imaginary part of the frequency 

We showed in the previous part that the conserved quanti- 
ties of a spherical photon orbit, {S',^,L^), correspond sim- 
ply to the real parts of the quasinormal-mode parameters, 
{cOR,Afij^,m), which are the leading-order quantities of a quasi- 
normal mode. Here, we will show that the behavior of the ho- 
moclinic orbits — namely, how the orbits propagate away from 
the spherical orbit, and how they move between 9± — reveals 
the spatiotemporal variation of the wave (i.e, the decay rate 
and the shape of its wave function in space). In Fig. 10, we 
plot the trajectory of a particular series of homoclinic orbits 
on the r-9 plane, to which we will refer at several points in 
the discussion below. 

With the appropriate values of ((?, ^,L^), the function u in 
Eq. (3.15) solves the wave equation to leading order and satis- 
fies the required boundary conditions. To recover the decaying 
behavior of quasinormal modes, however, we make correc- 
tions to the amplitude of the wave, which appear at next-to- 
leading order in the geometric-optics approximation. Because 
of symmetry, there should not be any correction to the ampli- 
tude in the direction, and the correction in the t direction 
should be a simple decay; therefore, we write 

M = Aexp(/5) = e-'''Arir)Ae{9)e-''^'e'^''^e^'^0e^'^' . 

v ' 

A{t.r.e) 

(3.24) 

This general expression contains four possible directions that 
the wave could be propagating: the ±0 direction and the ±r 
direction (depending on the signs in front of 5,- and Sq)- Be- 
cause the boundary conditions require that the waves propa- 
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-2 
log[(r-ro)/M] 

Figure 10: Schematic plot of trajectories in tiie r-d plane of lio- 
moclinic orbits outside of the peak of the potential (specifically for 
a black hole with spin a/M = 0.7 and a photon orbit with radius 
rg/M = 2.584). The two horizontal grid lines mark the turning 
points, = G±; between these turning points, there are two homo- 
clinic orbits passing through every point, while at turning points only 
one orbit passes through. Vertical grid lines indicate when the value 
of parameter A has changed along the orbit by (an arbitrarily chosen 
value) AA = 0.046M. Near the spherical photon orbit, each homo- 
clinic orbit undergoes an infinite number of periodic oscillations in 
6 while r — ro is growing exponentially as a function of X . 



gate towards +°o for r > ro and r^, — oo for r < ro, 
the sign in front of S,- should be positive for r > ro and neg- 
ative for r < ro- For motion, however, we insist that both 
directions (signs) be present, because a quasinormal mode is 
a standing wave in the 9 direction. Focusing on r > ro, we 
write 

= M++M_. (3.25) 

We will next require that both u+ and m_ satisfy the wave 
equation to next-to-leading order, separately. By explicitly 
computing Eq. (3.4) (or ^ const) in the Kerr spacetime, 

we find the amplitude satisfies the relation 



,^/logA 



d,{K{r)dX) 



1 



sin0 



(3.26) 

Here C, is an affine parameter along the geodesic specified by 
(ffji^i^j.). If we use the parameter X defined by d/dX = 
'Ld/dC, then we can separate the left-hand side of the equation 
as 



^d\ogA d , , , , d , , dt 



(3.27) 



Because the right-hand side of Eq. (3.13a) for dt/dX, sep- 
arates into a piece that depends only upon r and one that de- 
pends only upon 0, we will write Eq. (3.13a) schematically 



dt 
dX 



t + t. 



(3.28) 



where i is only a function of r and i is only a function of 9. 
Unlike in Eq. (3.13a), we will require that i average to zero 
when integrating over X for half a period of motion in the 6 
direction (i.e., from 0_ to 9+). We can ensure this condition 
is satisfied by subtracting an appropriate constant from i and 
adding it to i. Combining Eqs. (3.26)-(3.28) and performing 
a separation of variables, we obtain 



r—dlogA,- - 

iM — ^ -yt = — 

dr ' 4 

r-dlogA^ 



1 

2 sin 



(\/©sin0)', 



(3.29a) 
(3.29b) 



where a prime denotes a derivative with respect to r for func- 
tions of r only, and a derivative with respect to Q for functions 
of Q only (whether it is a or r derivative should be clear 
from the context). While it might at first seem possible to add 
a constant to the definition of i, and subtract it from i and still 
have both u+ and m_ satisfy the next-to-leading order geomet- 
ric optics, because we have akeady chosen to have i average 
to zero. 



(3.30) 



this separation is the only way to guarantee that |Ag | match 
each other at both ends. We will discuss the angular wave 
function in greater detail in the next part of this section. 

Let us now turn to the radial equation, from which we will 
be able to compute the decay rate. Close to ro, we can expand 
[%{r) to leading order as 



Substituting this result into Eq. (3.29a), we find 



dlogAr 
dr 



r-ro 




(3.31) 



(3.32) 



where we used the notation M'o = M"{ro). For A,, to be a func- 
tion that scales as A,- ^ (r — ro)" around ro for some integer 
n (namely it scales like a well-behaved function), we need to 
have 



= (n + 1 /2) lim 



dr jdX 



r^ro r — ro {dt/dX)g 



(3.33) 



To aiTive at the second line, we used Eq. (3.31), the fact that 
dr/dX = \fM, and that i is the part of dtjdX that does not 
vanish when averaging over one cycle of motion in the di- 
rection; the limit in the expression comes from the fact that the 
approximation in Eq. (3.31) becomes more accurate as r — > ro. 
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The physical interpretation of the rate that multiplies (n + 
1 /2) is somewhat subtle. Because the 6 motion is inde- 
pendent from r motion, a bundle of geodesies at the same r 
slightly larger than ro, but at different locations in 9, will re- 
turn to their respective initial values of 9 with a slightly in- 
creased value of r after one period of motion in the 9 direc- 
tion. The area of this bundle increases in the process, and by 
Eq. 3.8, the amplitude of the wave must decay; the rate of de- 
cay is governed by the quantity that multiplies (n + 1 /2) in 
Eq. (3.33). 

In addition, as shown in Fig. 10, the homoclinic orbits do 
pass through an infinite number of such oscillations near ro, 
because the radial motion is indefinitely slower than the 9 mo- 
tion as r approaches ro. It is clear from Fig. 10 that 



1 Ar Alog(r— ro) 



AA 



(3.34) 



approaches a constant as r — > ro. By multiplying the above 
equation by the constant value of (AA)/ (Af ) over one orbit of 
motion in the 9 direction. 



1 Ar 
r — ro Af 



Alog(r-ro) _ 

A^ = 



(3.35) 



also approaches a constant. This is usually defined as the Lya- 
punov exponent of one-dimensional motion; here, however, 
we emphasize that it is defined only after averaging over entire 
cycle of 9 motion. By comparing Eq. (3.35) with the second 
line of Eq. (3.33), and bearing in mind that the Lyapunov ex- 
ponent is defined after averaging over one period of 9 motion, 
one can write Eq. (3.33) as 



r= (n+2)7L- 



(3.36) 



To put Eq. (3.33) into a form that relates more clearly to 
Eq. (2.39), we use the conditions on the phase function. 



= 0, 



55 



= 0, 



(3.37) 



which hold for any point on the trajectory of the particle. We 
will apply this condition to two points on the particle's trajec- 
tory: one at (f , r, , ) and the second at (f + A? , r + Ar, , + 
A(^), where A/ is chosen such that the particle completes a cy- 
cle in 9 in this time (and it moves to a new location shifted 
Ar and A(/)). Substituting in the explicit expressions for the 
principal function in Eqs. (3.10) and (3. 1 la), we find 



d_ 
d 



A(r') 



= Af 



0, 



(3.38a) 
(3.38b) 



where we have defined 

^Se ^ 2 ^@{9')d9' = I ^@{9')d9' . (3.39) 



Because the change Ar is infinitesimal for r near ro, the in- 
tegrand is roughly constant, and the r-dependent part of the 



integral becomes the product of the integrand with Ar. Then, 
one can use Eq. (3.31) to write Eqs. (3.38a) and (3.38b) as 



1 



dM Ar ^^Sa 



yZ^Ao^^r-ro 

1 Ar <9A5, 



= Af, 



yZ^Ao^^r-ro 



(3.40a) 



(3.40b) 



Now, we also note that for a given fixed = m, the angular 
Bohr-Sommerfeld condition in Eq. (3.22) makes =S a function 
of S through the condition that ASg = (L— |m|);r. Because 
AS'e is a function of its total derivative with respect to <S 
must vanish. 



dASe dASe 



~dS 



= 0. 



(3.41) 



BS 



Therefore, when we multiply Eq. (3.40b) by {d£2 /dS)-s,s, and 
add it to Eq. (3.40a), we obtain the condition that 



1 



V^Ao 







(dB\ ' 


Ar 






KdS /bs. 





= At. (3.42) 



Combining this fact with the definition of the Lyapunov expo- 
nent in Eq. (3.35) and Eq. (3.36), we find that 



7 = 



V2^Ao 





5^ 


fd^\ 








KdS 





(3.43) 



where we recall that the quantities should be evaluated at ro. 
Equation (3.43) is equivalent to Eq. (2.39). Note, however, 
that in Eq. (3.43) we explicitly highlight the dependence of ^ 
on S through the term {d^ /dS')^^. There is an analogous 
term in Eq. (2.39) from the dependence of A/,„ on (O in the ex- 
pression for the potential V, which we must take into account 
when computing dV /dco; however, we did not write it out 
explicitly in Eq. (2.39). 

Summarizing the physical interpretation of the results in 
this subsection, we note that the Lyapunov exponent Yl is the 
rate at which the cross-sectional area of a bundle of homo- 
clinic rays expand, when averaged over one period of motion 
in the 9 direction in the vicinity of ro. The spatial Killing 
symmetry along (p means the extension of the ray bundle re- 
mains the same along that direction. This, therefore, allows us 
to write 

- e^'-' . (3.44) 

Correspondingly, the A\fs^ ~ const law requires that 

A^e-^^'l^, (3.45) 

which agrees with the decay rate of the least-damped QNM. 
The higher decay rates given by Eq. (3.33) come from an ef- 
fect related to the intrinsic expansion of the area of a phase 
front. More specifically, if the amplitude is akeady nonuni- 
form at points with different r — ro (but same 0), then shifting 
the spatial locations of the nonuniform distribution gives the 
appearance of additional decay. 
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3. Next-to-leading order: angular amplitude corrections and the 
imaginary part of Carter's constant 

Having found a relation in Eq. (3.29a) between the imagi- 
nary part of the energy, O)/, and the rate of divergence of rays, 
we now turn to Eq. (3.29b) to understand the geometric mean- 
ing of the complex part of A/„,. We recall from Sec. IIIC 1 
that i2 — Af^ — m^, at leading order, for a real Carter constant 
^. Because A/m becomes complex at next to leading order 
(and because m remains unchanged), if the correspondence 
^ = ^im — holds for a complex A/„,, then the Carter con- 
stant should also be complex, and its imaginary part should be 
equivalent to A|,j,. In this part, we argue that this relationship 
holds. 

By integrating Eq. (3.29b), we find that 



A± = 



1 



sin0V0 



exp 



± 







dO' 



(3.46) 



To interpret this equation, we will assume that the orbit is suf- 
ficiently close to ro that the change in r over the course of 
a period of motion in 9 is negligible. Under this assumption 
(and with the fact that dX=d9 / \/0) we can write the integral 
in the exponent in Eq. (3.46) as 



dB' 



=7 



- V0 
[f(0)-f(0-)]- 



Af 



[A(0)-A(0_)] 



(3.47) 



where Af/AA is the average of dt/dX over a cycle of Q mo- 
tion. We obtain this expression by using the fact that dt /dX 
is equivalent to i plus a constant when r (and hence f) does 
not change. Because i has zero average (by definition) over a 
period of Q motion, then when written in the form above, the 
constant must be {At) / (AA). We can write this average rate of 
change in a useful form by noting that, from Eq. (3.13a) and 
Eqs. (3.11b), 



dt 
c/A 



1 d^. 
lAdS 



-flVcos^e. 



(3.48) 



Averaging this expression over a cycle of motion, noting 
that the first term on the right-hand side is independent of 0, 
and using Eq. (3.11b) gives 



A«_ 
AA 



1 dS^, 
lAdS 

1 dS^ 
TABS 

1 dS^, 
lAdS 



e+ C0S2 



-dQ 



de 

7i 



dASe/dS' 
2dASe/d^ 
1 

2[dl,^, 



(3.49) 



In the last step we have used the Bohr-Sommerfeld condition 
(3.41). Also according to Eq. (3.13a) and Eq. (3.13c), we can 



find 



f(0)-r(0_) 



A(0)-A(0_) 



_d_ 



1 



2A dS 
d 



(A(0)-A(0_)), 



(3.50a) 
(3.50b) 



where to derive these two equations, we can again use the fact 
that dX =dd / \/® and the definition of 0; for the first we also 
make use of Eq. (3.48). 

Finally, we insert Eqs. (3.50a), (3.50b), and (3.49) into 
Eq. (3.47) to find 



- v© 



de' = (-/7) 



" d 




) -1 









[iSe{e)] . 



(3.51) 

Substituting Eq. (3.51) into the solution for Ag in Eq. (3.46) 
gives that 



Aj = 



exp < 



(3.52) 



The phase in this equation, however, is precisely the correc- 
tion to the leading-order expression for the phase e'^^ef^) jf 
allow S" and i? to be complex, where their imaginary parts are 
given by 



Iras' = -y=-(Oi, Im^ = 



Is 



(-7). (3.53) 



BS 



Through next-to-leading order, therefore, the portion of the 
wave is given by 



-AnC 



(3.54) 



where S and ^ used in Sg are complex. 

In the geometric -optics approximation, therefore, we have 
shown that we can account for the amplitude corrections to 
the wave by allowing the conserved quantities, S and to 
be complex [with their imaginary parts given in Eq. (3.53)]. 
Furthermore, through the geometric-optics correspondence, 
and the definition of A^^^^ in Eq. (2.27) we can confirm that 
A( = therefore, the relationship 



^^Ai,„-m^ 



(3.55) 



is true for a complex i? and A/,„. 

In closing, we note that at the same 0, the magnitude of the 
two components of the wave in Eq. (3.52) are not equal. More 
specifically, the integral involving f makes A+ have a larger 
amplitude at < n/2 and a smaller amplitude at > n/2; 
A^ has the opposite profile. Therefore, the net wave function 
remains symmetric about — n/2. 
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IV. FEATURES OF THE SPECTRA OF KERR BLACK 
HOLES 

In this section, we will use the WKB formula and the 
geometric -optics correspondence in the first two sections of 
this paper to explain several aspects of the quasinormal-mode 
spectrum of Kerr black holes. Specifically, we will explain 
the absence of damping for a significant fraction of modes of 
extremal Kerr holes. We will also decompose the frequency 
into orbital and precessional parts and explain a degeneracy in 
the spectra of Kerr holes in terms of a rational relation of these 
frequencies when the corresponding photon orbits close. 



A. Spherical photon orbits and extremal Kerr black holes 

We will first review the properties of spherical photon or- 
bits. These orbits can be found by setting ^(r) — SfJ{r) ~ 
(see, e.g., [34]), and their conserved quantities are fixed by the 
radius of the orbit r and the spin of the black hole a to be 



L-JS = 



{r' - 6Mr2 + 9M^r - 4a^M) 
fl2(r-M)2 



3Mr^ 



a{r~M) 



(4.1a) 
(4.1b) 



We will next discuss additional features of these orbits. 

For a given spin parameter a, there is a unique spherical 
photon orbit with parameters {<S,L^,£2) for any radius be- 
tween the outermost and innermost photon orbits (the ret- 
rograde and prograde equatorial photon orbits). Their radii 
(which we denote ri for retrograde and r2 for prograde orbits) 
are given by 



r2 =2M 



,2 f \a\ 

1 + cos I - arccos I ^ ^ 

,2 f\a\ 
1 + cos I - arccos I — 



(4.2a) 
(4.2b) 



At each ri < r < ri, the spherical orbit's inclination angle 
reaches a maximum and minimum of Q± (at which = 0). 
These angles are given by 

cos^0± = 



2^MA(2r3 - 3Mr2 +Mfl2) - {r' - m'^r + 2Ma^) 



(4.3) 



which are equivalent to the turning points of the integral (2.14) 
(and, therefore, we use the same symbols for these angles). 

Using the geometric-optics correspondence between 
((fjL^jcS) and [Hn^pL^af^^), we see that equatorial orbits 
at r\ and r2 corresponds to modes with /i = — 1 and +1, 
respectively, or modes with m = ±/ and I ^ 1 (strictly 
speaking, though, /i = m/{l + 1/2) never precisely reaches 
±1). We can also relate rp, the real root of Eq. (2.37), to 
the polar orbit and modes with m — Q. For orbits between 




Figure 11: The values of r and cos0+ of spherical spherical orbits, 
for a/M = (black, solid vertical line), 0.5 [red (light gray) dashed 
curve], 0.9 [blue (dark gray) dashed curve] and 0.99999 (black, solid 
curve). Note that for a = 0, all such orbits have r = 3M, while for 
a = M, a significant fraction reside at r = M. 



the equatorial and polar ones, we can use Eqs. (2.33a) and 
(2.33b) to obtain a ji between —1 and +1. Then, only those 
modes that can be written as m/{l + 1/2) with the allowed 
integer values of / and m correspond to a QNM (though the 
photon orbits that correspond to QNMs are a dense subset of 
all photon orbits). 

Note in Fig. 1 1 that for a ^ M, a significant fraction of 
spherical photon orbits of different inclination angles all have 
nearly the same radius, r w M. Through the geometric -optics 
correspondence, a large fraction of modes (a finite range of 
values of jj.) relate to this set of modes with M. In Fig. 12, 
we explicitly show the relation between modes character- 
ized by pL and their corresponding spherical-photon-orbit radii 
(normalized by the horizon radius) for several values of a/M 
slightly less than unity. The radius exhibits an interesting tran- 
sition between two kinds of behaviors; for > /i* w 0.744, 
the value of r is very close to M (the horizon radius for an 
extremal Kerr black hole), and for jj. < jj.^: the radii increase 
linearly. The orbits with ji > have a range of inclination 
angles. Their sin0± span from 0.731 (at /i*, the most inclined 
orbit) to 1 (at /i = 1, the prograde equatorial orbit). 

For the extremal black holes, therefore, a nonzero fraction 
of corotating spherical photon orbits appear to coincide with 
the horizon in the Boyer-Lindquist coordinate system. Al- 
though the proper distance between these orbits will not van- 
ish (see [36]), this does not seem to be a coordinate effect, be- 
cause there is a definite physical change of the modes for these 
values of jU > /i*. By comparing with Fig. 12 with Fig. 5, we 
see that these orbits also have i2/ ss 0. A vanishing imaginary 
part of the frequency corresponds to a vanishing of the radial 
Lyapunov exponent for this entire nonzero region of spher- 
ical photon orbits. This, therefore, would lead to a curious 
effect for a highly spinning black hole: for perturbations with 
fl > fl*, modes do not move away from or into the horizon 
very quickly. If we were to solve an initial-data problem con- 
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Figure 12: Radii of corotating spherical photon orbits as a function of 
/i, for a/M = 0.9 (black solid line), 0.99 (red dashed curve), 0.9999 
(blue dotted line). For extremal Kerr black holes, a nonzero fraction 
of all spherical photon orbits are on the horizon. 



taining these modes, we would find that they live for a long 
time; moreover, because these are nearly equatorial modes, 
we would see that the final, long-lived perturbations escape 
in the equatorial direction. This would imply that if we were 
to drive the black hole appropriately with equatorial incoming 
radiation, we would excite these nearly lossless modes with 
high amplitude (like optically exciting a resonant cavity with 
high finesse). 



B. A mode's orbital and precessional frequencies 

In this part, we will define two frequencies associated with 
individual spherical photon orbits (the orbital and precessional 
frequencies) and understand their connection to the real part of 
the QNM frequency. We begin by noting that because spher- 
ical photon orbits have only two independent degrees of free- 
dom describing their motion [see, e.g., Eq. (4.1)], the orbit 
can be characterized by two frequencies. The first is the 0- 
frequency, Q.q, the frequency at which the particle oscillates 
below and above the equatorial plane. During each 0-cycle, 
which takes an amount of time given by Tg = 2n/Q,g, the par- 
ticle also moves in the azimuthal (or (p) direction. If this angle 
is 2n for a corotating orbit (m > 0) or —2k for a counterrotat- 
ing orbit (m < 0), then there is no precession (and these sim- 
ple, closed orbits have effectively one frequency describing 
their motion, as the spherical photon orbits of a Schwarzschild 
black hole do). The difference between the A0 and ±.2n (its 
precession-free value) we will denote as the precession angle, 



A^prec = - 2;rsgnm , 



(4.4) 



where sgnm is the sign of m. We can also associate the rate of 
change of ^prec with a frequency. 



Hprec = A0prec/7e = A(j)prec^e / C^Tt) 



(4.5) 



Both Tg and A^piec can be computed from geodesic motion 
[see the formulas for Q.g and A^prec in Eq. (4.8)]. 



It is possible to perform split of the real part of the QNM 
into two analogous frequencies. To derive this split, start from 
a single ray, along which the phase of the wave must be con- 
stant. Also suppose that the ray originates from 0_ and ends 
at 0+ after traveling only one-half of a cycle of motion in the 
9 direction. During this time, the statement that the phase is 
unchanged is that 

= -(ORTg/2+{L-\m\)n + niA^/2. (4.6) 

Using (half of) Eq. (4.4), the real part of the frequency is 

(Or = LQ.g (m/L) + milprec (m/L) . (4.7) 

Note that Q.g and i2piec both depend on m/L. 

More explicitly, given the orbital parameters (S', ^,L^), the 
quantities Tg and A(j) can be obtained by computing 



(4.8a) 



Ad = - — 



1 f de 

2AdL^J^' 



(4.8b) 



(expressions that hold for any spherical photon orbit — not 
simply orbits that satisfy the Bohr-Sommerfeld condition) and 
the two frequencies are given by 



Q.g ^2n 



1 dM f de 

2aJIj 7i 



-1 



A(j) 



Oprec = £2e - (sgnL,)£2e . 

27Z 



(4.9a) 
(4.9b) 



These can be expressed in terms of {(S',J2,L,) using elliptic 
integrals (as was done in [37]), but we will not carry this out 
explicitly. 

For very slowly spinning black holes, a short calculation 
shows that 



(4.10a) 
(4.10b) 



where ro is the circular-photon-orbit radius for a 
Schwarzschild black hole, ro = 3M, and S — aM. The 
expression for Hg is the Keplerian frequency of the spherical 
photon orbit, and iiprec — ^S/tq is the Lense-Thirring pre- 
cessional frequency. In the slow-rotation limit, therefore, our 
formula recovers Ferrari and Mashhoon's result [9]. 

For any value of a, we can normalize Eq. (4.7) by L, and 
write 




Q.R{a,iJ.) ^Q.g{a,iJ.)+iJ.Q.precia,IJ.) 



(4.11) 



In Figs. 13 and 14, we explore the two frequencies in the de- 
composition of Q.R by separately plotting Q.g and iiprec, for 
different values of a. 
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Figure 13: Orbital frequency, Q.0, plotted against fl, for a/M = 0.3 
[red (light gray) solid curve], 0.7 [blue (dark gray) solid curve], 0.9 
(purple dashed line), and 1 (black dotted line). The orbital frequency 
vanishes for a significant range of fl for extremal black holes. 




Figure 14: Precessional frequency, fl0, versus /x plotted identically 
to those curves in Fig. 13 representing the same black-hole spins. 
The precessional frequency approaches the horizon frequency, fl//, 
for a range of values of /i for extremal black holes. 



For small values of a/M, Q.g and Hpiec are consistent with 
the constant values predicted by Eqs. (4. 10a) and (4. 10b). For 
larger values of a/M, Q.g does not vary much as a function 
of jj. until a ^ Q.IM; for spins greater than this value, it is 
only for larger values of jj. that Q.g changes significantly by 
decreasing from the equivalent values for a = 0. Finally, as 
a M, Q.g vanishes for > /x* w 0.744. The precessional 
frequency, £2prec, on the other hand, monotonically increases 
as a function of jj, even for small values of a/M; as a — > M, 
£2piec grows to Q.H at ^ /i*, and stays there for all values of 
/i > /i*. For a^ M and jj. > there is one additional feature 
worth noting: because Q.g ^ and Q.^ ^ Q.^, this gives rise 
to the interpretation of the mode as a ray that sticks on the 
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Figure 15: A diagram showing the spin parameters, a, and the ratios 
of the multipolar indexes in/L, at which the orbital and precessional 
frequencies have a ratio of p/q. Although we only perform our nu- 
merical calculations at a discrete set of m/L values (shown by the 
dots), in the eikonal limit, each set of points for a given ratio of p/q 
approaches a continuous curve. 



horizon and corotates with the horizon at its angular frequency 
of Q.H', moreover, there seems to be no restoring force along 
the 9 direction. 



C. Degenerate quasinormal modes and closed spherical 
photon orbits 

Finally, in this section, we interpret the degeneracy of QNM 
frequencies (of which Fig. 1 was an example). Recall that 
in that figure, for a/M w 0.7, we found pairs of modes such 
as (2,2) and (3,-2), (3,2) and (4,-2), (4,2) and (5,-2), 
and so on, all have approximately the same frequency. For 
another, lower spin a/M « 0.4, pairs like (3,3) and (4,-3), 
(4, 3) and (5,-3), et cetera, have approximately the same fre- 
quency. 

The approximate degeneracy exists because the ratio be- 
tween Q.g and iipiec can be rational, and the photon orbits 
close. If for a certain mode of a black hole with spin a, with 
m and L, and for integers p and q, 

(in \ / m\ 

a,j-)^P^'ec[a,j-)^ (4.12) 

this means that there exists a closed spherical photon orbit that 
satisfies the conditions necessary to coiTespond to a QNM. 
Equation (4.12) implies that 
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If Clg and Hpiec do not change much from jj. = m/L to ji' = 
[m — kp) /{L + kq) (either because spin is small — and there- 
fore He and Hpiec depend weakly on /i — or because kq 
and m ^ kp), then 
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(4.14) 
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Figure 16: For black holes with spins a/M = 0.768, 0.612, and 0.502, the spherical photon orbits with Ogrb = 2a)prec, on the left, DJorb = Sfflprec 
in the center, and atath = 4tt)prec on the right, respectively. These orbits correspond to quasinormal modes in the eikonal limit with m/L = 0.5. 
The top figures show the photon orbit, the red, solid curve, on its photon sphere (represented by a transparent sphere). The dashed black line is 
the equatorial (0 = 7t/l) plane, which was inserted for reference. The bottom figures are the same photon orbits, but plotted in the <p-G plane, 
instead. 



Because Clj depends similarly on /z, under the same condi- 
tions, 
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(4.15) 



therefore, the modes are degenerate. It is also clear from 
Eq. (4.12) that the degeneracy happens at the same time that 
the corresponding orbit is closed. The three series mentioned 
at the beginning of the paper correspond to p/q = 4,6, and 8, 
respectively (for k — 1). 



1. Slowly spinning black holes 

For a/M < 1, when Eqs. (4.10a) and (4.10b) apply, the 
condition for degenerate modes becomes 



qV27 



M 



< 1 



(4.16) 



(a statement that holds independent of /x). This implies that 
orbits of all inclinations close for these spins. 

For these specific spins, the QNM spectrum is completely 
degenerate, by which we mean that all modes have the same 
decay rate, and all real parts of the frequencies are inte- 
ger multiples of only one frequency (similar to those of a 
Schwarzschild black hole). Using this approximate formula 
to find a for the three instances of degeneracy in Fig. 1, we 
find 



2. Generic black holes 

For a generic spin parameter a, we will explain degenera- 
cies that exist around a mode with L ^ 1 and \m\ ^ 1. If the 
condition in Eq. (4. 12) holds for p,q<^ min(L, |m|), then there 
is a range of \k\ ^ min(L, \m\)/ max{p,q) in which there is a 
degeneracy between all {L + kq,m — kp) and (L,m). These 
modes must be those close to the mode of indices {L,m), be- 
cause, strictly speaking, it is only the orbit coiTesponding to 
m/L which is precisely closed. 

To find this degeneracy, we will search for spin parameters 
a for which Eq. (4. 12) holds for any set of indexes (L, m) and 
integers {p,q) that satisfy L,\m\ ^ p,q (we generally either 
find one or zero solutions). To visualize this degeneracy, for 
each pair {p,q), we will mark all possible pairs of {m/L, a) 
in a 2D plot; the values of the spins are sufficiently dense for 
each value of m/L that they form a smooth curve when plotted 
against m/L. Some of these curves are shown in Fig. 15. Be- 
cause for a fixed p/q the degenerate spins for a/M < 0.3 are 
nearly independent of m/L, Eq. (4.16) should be an accurate 
prediction for spins less that that value. As a concrete illus- 
tration of the orbits corresponding to these degenerate modes, 
we plot closed orbits for m/L — 0.5 and for a/M « 0.5, 0.61, 
and 0.77 in Fig. 16. The values of the spins agree quite well 
with those predicted in Fig. 15. 



V. CONCLUSIONS AND DISCUSSION 



fl4/iw0.65M, fl6/iw0.43M, fl8/i«0.32M. (4.17) 
These are not very far away from spins we found empirically. 



In this paper, we extended the results of several earlier 
works [9, 27, 31, 32] to compute the quasinormal-mode fre- 
quencies and wave functions of a Kerr black hole of arbitrary 
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astrophysical spins, in the eikonal limit (I ^ 1). We focused 
on developing a greater intuitive understanding of their behav- 
ior, but, in the process, we calculated expressions for large- 
/ quasinormal-mode frequencies that are reasonably accurate 
even at low I. Specifically, we applied a WKB analysis to the 
system of equations defined by the angular and radial Teukol- 
sky equations. Using a Bohr-Sommerfeld condition for the 
angular equation, we related the angular separation constant 
to the frequency; when we expanded the constraint to lead- 
ing order in aco/l, we found an equally accurate and alge- 
braically simpler relation between the frequency and angular 
eigenvalue. We then used a well-known WKB analysis on the 
radial Teukolsky equation to obtain expressions for the QNM 
frequencies and the angular separation constants. The accu- 
racy of the approximate expressions for the QNM frequency 
are observed to be of order 0{L^^) even though we had only 
expected a (9(L^^) convergence for the imaginary part. 

Next, we reviewed that a massless scalar wave in the 
leading-order, geometric -optics approximation obeys the 
Hamilton-Jacobi equations, which are very similar to the 
Teukolsky equations when Z ^ 1. By identifying terms in 
the Hamilton-Jacobi equations and Teukolsky equations, we 
related the conserved quantities of the Hamilton-Jacobi equa- 
tions to the eigenvalues of the separated Teukolsky equations. 
Specifically, we confirmed that the energy, angular momen- 
tum in the z direction, and Carter constant in the Hamilton- 
Jacobi equations correspond to the real frequency, the index 
m, and the angular eigenvalue minus in the Teukolsky 
equations, respectively. Furthermore, we found that the con- 
ditions that define a quasinormal mode in the WKB approxi- 
mation are equivalent to the conditions in the geometric-optics 
approximation that determine a spherical photon orbit that sat- 
isfies an identical Bohr-Sommerfeld quantization condition. 

By analyzing the next-to-leading-order, geometric -optics 
approximation, we showed that the corrections to the ampli- 
tude of the scalar wave correspond to the imaginary parts of 
the WKB quantities. Specifically, we saw that the imaginary 
part of the frequency is equal to a positive half-integer times 
the Lyapunov exponent averaged over a period of motion in 
the Q direction. The imaginary part of the angular eigenvalue 
is equal to the imaginary part of the Carter constant, which is, 
in turn, related to an amplitude correction to geometric -optics 
approximation to the angular function for 0. 

We then applied these results to study properties of the 
QNM spectra of Kerr black holes. We observed that for ex- 
tremal Kerr black holes a significant fraction of the QNMs 
have nearly zero imaginary part (vanishing damping) and their 
corresponding spherical photon orbits are stuck on the hori- 
zon (in Boyer-Lindquist coordinates). We plan to study this 
unusual feature of extremal Kerr black holes in future work. 
In addition, we showed that for Kerr black holes of any spin, 
the modes' frequencies (in the eikonal limit) are a linear com- 
bination of the orbital and precession frequencies of the cor- 
responding spherical photon orbits. This allows us to study 
an intriguing feature of the QNM spectrum: namely, when the 
orbital and precession frequencies are rationally related — i.e, 
when the spherical photon orbits are closed — then the corre- 
sponding quasinormal-mode frequencies are also degenerate. 



We hope that the approximate expressions for the 
quasinormal-mode frequencies in this paper will prove help- 
ful for understanding wave propagation in the Kerr spacetime. 
This not unreasonable to suppose, because Dolan and Ottewill 
have shown in [26, 28] that to calculate the Green's func- 
tion analytically in the Schwarzschild spacetime, one needs to 
know analytical expressions for the frequency of the quasinor- 
mal modes (specifically, this comes from the fact that the fre- 
quencies of the quasinormal mode are the poles of the Green's 
function in the frequency domain). We, therefore, think that 
our approximate formulas could assist with the calculation of 
the Green's function in the Kerr spacetime, in future work. 
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Appendix A: The Taylor expanded Bohr-Sommerfeld condition 

The Bohr-Sommerfeld constraint (2.14) gives us a way to 
evaluate A/,„ in terms of /, m, and o; the error in this ap- 
proximation scales as 1 //. Because it is an integral equation, 
it is not particularly convenient to solve, and it is beneficial 
to have an approximate, but algebraic expression for the fre- 
quency of a QNM. With the benefit of hindsight, one can 
confirm through numerical calculations of exact QNM fre- 
quencies performed using Leaver's method that the parameter 
a(o/l is numerically a small number for all black hole spins. 
We can then expand the angular separation constant, A/„,, in 
a series in a(o/l as A/,„ = A^J^^ + 5A/,„, where A^^^^ satisfies the 
equation 

and at leading order, Q^j* , 0" = ± arcsin[m/ (/ + 1 /2)]. One can 
easily verify that the solution to this equation is the angular 
eigenvalue of a Schwarzschild black hole, A^j^^ — + 1/2)^ 
(note that we are assuming / ^ 1). Now we will compute 
the lowest-order perturbation in ad) /I, which turns out to be 
quadratic in this parameter [i.e., (ao)//)^] below: 

^A.„ + .^a.^cos^e ^^^^ ^^^^ 
•y(Z + l/2)2-m2/sin2 
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The integration limits 0+ , 0_ also can be expanded in a se- 
ries in flO)//, and the lowest-order terms of this series are given 
by 0^,0''; The perturbation in would result in some 

quartic corrections in a(o/l [i.e., (ao)//)^] when we evaluate 
the integrals of Eqs. (A2) and (Al), because the integrand is 
of order (aco/l)^ and the width of the correction in 0+, 0_ are 
also of order {a(o/l)^. As a result, we will not need it here. 
Evaluating the integral in Eq. (A2) is straightforward, and we 
find 



9 9 
a-(0 



1- 



(A3) 



Interestingly, the above expression is consistent with the ex- 
pansion of A/,„ for small ad) given in [39], even in the eikonal 



limit, where aco is large. The reason for this fortuitous agree- 
ment is again that for QNMs of Kerr black holes of any spin, 
ad) /I is small, and the expansion only involves even powers 
of this parameter, {a(o/l)^. Because the coefficients in the 
expansion of a(0 scale as 1 /l'' for even powers of (00)*^ and 
1 for odd powers of (aO))^, in the limit of large I, the two 
series actually are equivalent in the eikonal limit. In principle, 
one can also expand and solve Eq. (2.14) to higher orders in 
the parameter aco /I and compare with the expansion in a o in 
[39]; we expect that the two series should agree. This would 
be useful, because it would effectively let one use the small 
aco expansion in the eikonal limit, where the series would, os- 
tensibly, not be valid. 
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